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We propose a spin transport experiment to measure the low-energy excitations in insulators 
with spin degrees of freedom, with a focus on detecting ground states that lack magnetic 
order. A general formalism to compute the spin-current from a metal with a non-equilibrium 
distribution of spins to an insulator is developed. It is applied to insulating states with and 
without long range magnetic order, and salient features in the spin-conductance are noted. 
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I. INTRODUCTION 

Observation of fractionalized excitations in insulating spin-systems has been a long-sought goal 
in physics. Such quantum spin liquid states, if realized in nature, would be a new quantum phase of 
matter with exotic properties. Certain candidate materials have strong experimental evidence for 
exhibiting spin liquid ground states. For example, thermal conductivity experiments on insulating 
frustrated triangular lattice organic salts by M. Yamashita et al. indicate presence of mobile 
gapless excitations. Inelastic neutron scattering experiments on single crystals of Herbertsmithite, 
a kagome lattice spin-half Heisenberg antiferromagnet by Han et al. provide evidence for the 
presence of a continuum of fractionalized spinon excitations. Numerical studies on the triangular ’ 
and kagome 1 lattice Heisenberg models also indicate the possibility of spin liquid ground states in 
certain parameter-regimes. 

In spite of promising evidence for observation of spin liquids from several experiments 1,2,<v , the 
exact nature of experimentally realized ground states, and in particular, the presence of a spin-gap 
is still unclear. In this paper, we propose a transport experiment which can probe the mobile 
spin-carrying excitations of the system at low energies; these experiments are similar in spirit 
to those discussed recently by Takei et al. ~ and collaborators L ~ for materials with magnetic 
order. Recent advances in spintronics 1 ’ have made it possible to create a spin-accumulation at 
boundaries of metals via the spin Hall effect. We propose to use this non-equilibrium accumulation 
of spins to inject a spin-current into an insulating state with spin-degrees of freedom. The spin- 
current is a function of the spin-accumulation voltage in the metal. Therefore, by measuring 
the spin-current as a function of this voltage, and looking at thresholds and exponents, we can 
comment on the presence of spin-gaps and the low-energy dispersion of the fractionalized spin-half 
excitations. 

The rest of the paper is organized as follows. In section II, we describe the geometry of our 
setup, and develop a formalism to evaluate the spin-current injected into a magnetic insulator from 
a metal. In section III, we apply the formalism to evaluate the spin-current into an antiferromagnet 
with collinear Neel order. In section IV, we first analytically calculate for the spin-current into 
insulating states with no long range magnetic order, including both valence bond solid states and 
spin liquid states. Then we go beyond the analytical approximations, and numerically identify 
some broad features in the spin-conductance for a spin liquid ground state on the kagome lattice, 
which is a candidate state for Herbertsmithite 2,1 '. Details of relevant calculations are contained in 
the appendices. 
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II. FORMALISM TO EVALUATE SPIN-CURRENT 

A. Generation and detection of spin-current 


We begin with a brief discussion of the spin Hall effects, which we shall use to generate and 
detect spin-currents, and then describe the exact geometry of spin injector and detector we use. A 
charge current passed through a paramagnetic material can drive a transverse spin current in pres¬ 
ence of strong intrinsic spin-orbit coupling or skew-scattering by spin-orbit coupled disorder “ . 
The spin current impinging on the boundary is given by Js = t^OshJc, where Jc is the charge 
current density and Osh is the spin-Hall angle, and sets up a spin-accumulation at the boundary, 
that has been measured in experiments for both metals 22, and semiconductors 2 “ . The recip¬ 
rocal process, where injecting a spin current into a spin-orbit coupled paramagnetic material sets 
up a charge current (or voltage) transverse to the spin-current - the inverse spin Hall effect, has 
also been observed ,26,2 '. Furthermore, both processes have been used simultaneously to transmit 
electrical signals across a magnetic insulator . Theoretical predictions for the spin superfluid trans¬ 
port through a ferromagnetic lS and antiferromagnetic insulator sandwiched between two metallic 
reservoirs have been worked out in the linear response regime. Taking phenomenological Gilbert 
damping into account, the spin current density J r s pumped into the right reservoir as a function of 
the spin accumulation voltage V is given by 8,!) 


T r — 
J S — 


V 


n u 

9l 9r 


+#+<,„ 


a) 


where gj^ is the spin flip conductance at the left (right) interface, and g a quantifies the loss in 
spin current due to Gilbert damping. 


Spin accumulation 



(a) Spin accumulation via the spin Hall effect, 
and injection at the left interface 



(b) Spin-current detection via the inverse spin Hall 
effect in the right metallic reservoir 


FIG. 1: Geometry for generation and detection of spin-current 
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Let us consider an analogous geometry, where an insulating block with spin degrees of freedom 
is placed in between two metallic reservoirs, as shown in Fig. 1. A charge current in the left metallic 
reservoir, in presence of strong spin-orbit coupling, will create a non-equilibrium accumulation of 
spin at the metal-insulator boundary. We assume that there are no thermal gradients, and that the 
spin accumulation can be well modeled by different chemical potentials /x-j- and /x^. in the Fermi Dirac 
distribution at temperature T for the spin-up and spin-down electrons. The left metal reservoir 
will subsequently relax by sending a spin-current into the spin insulator. We assume negligible 
loss of spin-current inside the insulator, so that the spin-current sets up a spin-accumulation at 
the insulator-metal boundary on the right. If the metallic reservoir on the right was initially in 
thermal equilibrium at T, the accumulated spin density at the boundary will drive a charge current 
via the inverse spin Hall effect. This charge-current, or the associated voltage can be detected, and 
therefore we can find the spin-current by measuring charge currents (or voltages) in both metallic 
reservoirs. 


B. General expression for spin-current 

Let us choose x as the longitudinal direction which is normal to the interfaces, and z as the 
spin-quantization axis. We shall evaluate the spin-current crossing the left metal-insulator interface 
when V = /x-j- — /xj_ > 0. To make analytical progress, we assume a clean interface between the 
metal and the insulator, with translational invariance in the plane of the interface. The metallic 
reservoir is assumed to be a Fermi liquid with quadratic dispersion and Fermi energy ep, so that 
n CT (e) = ^ + 1^ with (setting h = 1). We shall always work in the regime where 

T, V <C €f, and henceforth set /x-j- = /x, so that /xj, = /x — V, to simplify notations. 

We assume that the electron spin S e in the metal interacts with the boundary spins of the 
insulator, located at interface lattice sites Xj. via a local spin-rotation symmetric local Hamiltonian 

H int = jJ2Se-S j 8(x e -X j ) ( 2 ) 

3 

Let the insulator have exact eigenstates (|n)}, then its initial state is described by the equilibrium 
density matrix e ^ n | n) (n\. For the metal, periodic boundary conditions in a large box of 
volume V = L X A± is assumed, where A± is the interface area. We now use Fermi’s golden rule 
to calculate the rate of scattering of a right-moving electron state |£q,t) to a left-moving electron 
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state \k 2 ,\.}. The matrix element for scattering to a final state | m) of the insulator is given by 


J 


(k 2 , m\ H int |*h, t; n ) = ^ E e ^' Xj ( m l S t ^ ’ definin § q = h~k 2 


(3) 


Defining u(ki, k 2 ) = ^ ^ as the energy transfer, the rate of scattering R is 


m,n 
7tJ 2 


r 


—/3E n 


‘ILIA i 


S-+ q±,u = 


(fci,t;n| H int \k 2 ,i;m) 
2ki ■ q - q 2 \ 


& ( -^n m € 


k2,i 


2m 


(4) 


where S _[-('LL; w) is the dynamic spin structure factor of the insulator at the interface, defined as 


S-+(q±,u,) = — E j df e i„t {Sr{t) s+( 0)) thermal 


(5) 


hi 


The spin-current crossing the boundary for this scattering event is If we have R such events 
per unit time, then the net spin-current crossing the boundary is just Summing over all initial 
electron and final states consistent with phase space constraints, the current I sp in,t due to up-spin 
electrons getting reflected to down-spin ones is 


/, 


spin 


7T u A _\_ 

4 m 


d d k i 


d d q 


Jk\ x >o (2ir) d J q x >k\ x ( 2 tt ) 


I n t(%) i 1 - n f(%-g)) Qx S-+ j ?X,W = 


2k ■ q — if 


2m 


( 6 ) 


At non-zero T, the reverse process where spin-down electrons get reflected to spin-up ones con¬ 
tribute analogously a spin-current I sp in,i given by 

nJ 2 A± f d d k\ f 


I, 


spin , 4 - 


4m Jk lx > o (2vr) d J qx>klx (2vr) 


d dq , s( ( A a (- 2k-q-q 2 

d n f(%) ( 1 - M e k,-a) ) 5 +- I w = 


"ki—q 


2m 


The net spin-current is therefore given by the difference of the two contributions listed above 


Run, — Isriiril T 


spin — -l spin, i~ •‘spin ,4- 


(7) 


( 8 ) 


C. Simplifications for certain physically relevant structure factors 

The expression for the spin-current can be considerably simplified once we note that at T —> 0, 
scattering is essentially restricted within an energy window of V. For oj < V, we assume that the 

dynamic structure factor S. t _ (q±,u>) assumes large values only for small This is physically 

relevant for several systems where excitations at large momenta typically have large energy cost. 
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FIG. 2: Allowed phase space for scattering of an electron with given initial momentum 


As Fig. 2 shows, if the system does not have excitations at oj < V for \q± \ > A, then scattering is 
restricted within a patch of dimensions ^ X A d ~\ v F being the Fermi velocity. 

To exploit this, we approximate the initial momentum k\ ~ k F h, and linearize the energy 
transfer u about the point of elastic scattering as follows 


q = 2 k F {h ■ x ) — 5q x x — q± 
u(ki,q) = v F [(n • x)Sq x - h ■ q±\ + 0{5ql, q\) 


(9) 


We also assume that the electronic density of states v(e F ) is approximately a constant near the 
Fermi surface for Sq x ,q± <^k F - Leaving the details of calculation to appendix A, these simplifi¬ 
cations lead to the following form of the spin-current for spin-up electrons flipping to spin-down 
ones. 


T _irJ 2 A ± v(e F ) f du d d l q x {V - u) c r _, ^ 

spin,t “ 4 J 2n (2Tr) d ~ 1 1-e-^ v - LJ >~ +[q±,U}) 

Analogous manipulations for the reverse process lead to 

T _ irJ 2 A ±F (e F ) f du d d ~ l q x (V + u) Q ( _ ^ 

lspin ’ 1 ~ 4 J 2n (27r) d - 1 efA v +“) - 1 + ~ [q±,U) 


( 10 ) 


( 11 ) 


These expressions make it transparent that as T —> 0, only up-spin electrons flipping to down-spin 
ones contribute the energy window (0,17). The reverse process is always exponentially suppressed 
as there must be an energy gain of at least V for a down-spin electron to flip to an up-spin one due 
to phase space constraints. The net spin-current is, as described in equation (8), the difference of 
the above two currents. 
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This formalism can be extended to cases where the quasiparticle excitation energy has minima 
at large transverse momenta {Q±} (with magnitude of a -1 where a is the microscopic lattice 
length-scale), provided the different Q± are well-separated from each other. This is typically true 
for systems with quasiparticle bands, as the momenta difference between the band minima are of 
the order of a~ 1 . For example, cubic lattice antiferromagnets with a 2 dimensional boundary have 
spin-wave excitations about the ordering wave-vector ^ (0,1,1). Referring the reader to 

appendix A again for the details, here we just state the main result. The effect of inelastic scattering 
about large transverse momenta Q± is to scale the spin-current by an overall 0( 1) angular factor 
fang(kp/Q±), so that equation (10) for I sp i n -f is now modified to 


/, 


spin ,t 


7T J 2 Aj_i/(e F ) 


Y fang{k F /QJ_) 
Q± 


duo d d 1 q F (V — uo) 

2ir (27r) rf_1 1 — e—P( v ~ u ) 


S- + tf ± ,u) (12) 


where the angular factor, coming from kinematical constraints, is given by 

( 


fang(k F /Q±) — 


' k 2 F (h-x) 2 +'2k F {Q_L-h)>Q 2 L 


dQ 
Sdl i 


i + 


k F (h ■ x) 


y ( k 2 F (n ■ x) 2 + 2k F (Q± ■ n) - Q\ 


1/2 


(13) 


In equation (13), Sd -1 is the sphere in M rf , and one can check that for Q± = 0 the angular factor 
reduces to unity, as desired. One can also check the limit Q_l k F , in which case scattering of 
the electron by q± « Q± is excluded by phase space constraints and f a ng{k F /Q± ) —^ 0. Equation 
(11) also undergo similar modifications, and putting these together we obtain our main result of 
this section 


Ispin — 


irJ 2 A± F (e F ) 


T. fang{k F /Ql_) 

Qx 


duo d d 


-i 


q± 


2n (27r 


\d— 1 


{V-u) c ^ A (V + co) 0 ^ A 

i _ e -p{v-u}) (i± ; w) e p(v+u) _ i s+- 


(14) 


We once again carefully note that this formalism for extension of the spin-current calculation to 
a set of different {Q_l} works only when the different points are well-isolated in the Brillouin zone 
of spin-carrying excitations of the insulator. Physically, this implies that the different momentum 
patches (to which the electron is scattered) do not overlap with each other. If they start to overlap, 
then we would count the same final electron state multiple times and over-estimate the spin-current. 


III. SPIN CURRENT FOR ORDERED ANTIFERROMAGNETS 

In this section, we apply the formalism developed in section II to calculate the spin-current from 
the metallic reservoir to an ordered collinear antiferromagnet, deep in the Neel phase. We assume 
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d = 3, so that a symmetry-broken state can occur at T > 0. The results can also be generalized 
to d = 2 at T = 0. In the following subsections, we illustrate evaluation of the current with the 
simplest scenario - a cubic lattice antiferromagnet with ordering wave vector Q AF = ^ (1,1,1), 
so that Q AF = ^ (0,1,1). We split our analysis into two subsections, corresponding to the Neel 
order pointing perpendicular and parallel to the spin-quantization axis in the metal, and add up 
the contributions due to elastic reflection from the static magnetic moments, and the inelastic 
reflection due to spin-wave excitations, to find the net spin-current. 


A. Neel order perpendicular to spin quantization axis in the metal 

1. Elastic contribution 


In order to contribute the elastic spin-flip scattering from the metal-antiferromagnet interface, 
we replace the fluctuating spin operators at the boundary by static moments, resembling the 
classical ground state. For Neel order along y, which is normal to the spin-quantization axis z in 
the metal reservoir, we can write the Hamiltonian as 


H int = J^2s e -Sj 6(x -Xj)^J S J2 s y S(x - Xf) (15) 

3 3 

We use Fermi’s golden rule again to find the rate of scattering of spin-flip scattering of electrons 
at the interface 


R = 2ir| <*a,4.| H,„t |*i, t) I y(% - %) = 


2s,. . , _ ^ j. 


5 ( % - %-q 


(16) 


4L2 

Following an analogous procedure of finding the spin current due to this scattering event, and 
summing over all initial and final states consistent with phase space restrictions, we arrive at the 
following expression for the elastic contribution I sp m in terms of f an g (k F /Q AF ) 


jel 

‘-spin, 


jel 

‘spin, 


lt = fang (iWQl ) 


AF\ kJ 2 A ± v(e F )V 


X = fang ( k F /Qj _ F ) 


8 1 - e-W 

AF\ nJ 2 A ± u{e F )V 


ePV - 1 


Ifpin = I spin,X - I spin,X = fang () 


AF\ Fj 2 A L v{e F ) 


V 


(17) 

(18) 
(19) 


Note that the elastic contribution to the current is proportional to the number of propagating 
modes at the Fermi surface, given by v(e F )V. So this contribution is similar to what one would 
obtain by using the Landauer formalism, as had been done for an analogous geometry by Takei et 
aV . 






10 


2. Inelastic contribution 


The inelastic contribution can be directly evaluated by application of equation (14), as the 
ordered antiferromagnet deep in the Neel phase has spin-wave excitations that have minimum 
energy about Q± = 0 and Q± = Q‘ 4F , which are well-separated in the insulator Brillouin zone. We 
work in the T —> 0 limit, which implies that the insulator is initially in its ground state. Therefore 
oj > 0 in the dynamic structure factors, and we can drop the contribution from ^ to the 
spin-current. 

We use the Holstein-Primakoff transformation to diagonalize the Hamiltonian and evaluate 

S _ (q±,u). Leaving the details to appendix B, the dynamic structure factor in the small and 

T —> 0 limit is given by (for cj > 0, setting a = 1) 


S-+ (q ± ,u) 


7T q± 
8\/2 


5 (uj 


v s q±) 


( 20 ) 


where v s is the speed of spin-waves in the antiferromagnet. We can plug this back into equation 
(14), and we obtain the inelastic contribution to be 


Jinel T->0 irj‘ 2 Aj_u(e F ) 
spin a 


[1 + fang ( k F /Q f F )] 


V 4 

384\Z2ttv | 


( 21 ) 


We now add up the contributions from equations (19) and (21) to find the net spin-current 
when the Neel order is perpendicular to the spin-quantization axis in the metal. 


T—>o irJ 2 Aj_ F (e F ) 


L spm 


fang {kp/Qf) V + [l + fang {k p/Q^)} 


V 4 


192 ^ 2 - 


■KVi 


( 22 ) 


B. Neel order parallel to spin quantization axis in the metal 

1. Elastic contribution 

For Neel order along z, which is normal to the spin-quantization axis z in the metal reservoir, 
we can write the Hamiltonian as 

H int = J^S e -Sj 5(x - Xj) ->■ JS S z s(x - Xj) (23) 

3 3 

In this case, the Hamiltonian H m i commutes with the z-component of the electron spin, and 
therefore cannot flip it. Therefore there is no elastic contribution to the spin-current. 
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2. Inelastic contribution 

For the inelastic contribution, we again use the T —> 0 limit of equation (14). The dynamic 
structure factor is evaluated in an analogous manner to the previous subsection III A 2, and is 
essentially identical to equation (20) barring a constant extra pre-factor. We find that the net spin 
current when the Neel vector is along the spin-quantization axis is given by 

I** T =° cjl T =° Ij2 ' 4 f (eF> [i + u, (wq± f )] (24) 

IV. SPIN CURRENT FOR SYSTEMS WITH NO MAGNETIC ORDER 

In this section, we shall apply the formalism from section II to evaluate the spin-current into 
states with no long range magnetic order. Some candidate phases for Mott insulators with unbroken 
spin-rotation symmetry are described by spin-half quasiparticles or spinons, coupled to an emergent 
gauge field. In the deconfined phase of the gauge field, the lattice symmetry is unbroken and the 
ground state is a spin liquid . The spinons can propagate as independent quasiparticles and carry 
a spin-current. In the confined phase, the ground state might spontaneously break translation 
symmetry of the lattice, resulting in a valence bond solid (VBS) state with short-range order. 
In this case, the low-lying excitations with non-zero spin are spin-triplets or triplons, which are 
gapped excitations that carry the spin current. 


A. VBS states with triplon excitations 


At low energies, the structure factor will be dominated by single triplon excitations. Let us 
assume that the triplon has a gap A^ and a quadratic dispersion, so the dynamic structure factor 
can be approximated by 

S-+ (g_L, to) « C 5 (u - A t - qgj 2 ) (25) 


Here we also assume that the prefactor C is independent of u> and q±. Now we again use the 
T —> 0 limit of equation (14) to compute the spin-current. For a d dimensional system with ad-1 
dimensional boundary, we find that the spin-current is given by 


L spin 


T A o {v _ A \d/2+i e(v _ A ) 

2(2ir) d d(d + 1) 1 1 1 


(26) 


As expected, there is a threshold at V = At, as energy conservation implies that no triplons can 
be excited when V is less than the triplon gap. Above the cutoff, the spin-current has a power law 
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behavior with voltage with an exponent that depends on the dimensionality d of the system. For 
instance, in d = 3, the exponent is |. 


B. Spin liquids with spinon excitations 


We first approach the problem analytically by using a low energy effective theory to calculate the 
two-spinon structure factor. We use a mean-held approach where the spinons are free quasiparticles 
in the system, and have negligible coupling to other excitations which do not carry spin (like visons, 
which are vortices of the emergent gauge held). For a given spinon dispersion eg, the free-spinon 
Green’s function in imaginary time is given by 

G s (k,iu n ) = — - (27) 

KO n - €g 

where oj n is a Matsubara frequency which is determined by bosonic or fermionic statistics of the 
spinons. We can calculate the structure factor from the dynamic susceptibility x _h given by 


X — ^ ' G s ( k , i£ln)G s (k -\- q± , i£l n i^n) 

k,itt n 

f d 2 k A -n B (eg)-n B (eg +o - i )' 


T—s>0 


(271r) 2 
d 2 k 


-iu n + eg + e k+?± 


k+q± ] (for bosonic spinons) 


(2vr) 2 (-iwn + eg + eg^J 


(28) 


which, in turn, leads to the following result for the zero-temperature limit of the dynamic structure 
factor 


S-+(q±,u) = 1 _ e -p u Im [x-+(gl, iu n u + irjj] 

T ~°^ >0 hm lm[x- + (q±,iu n -+ u + ir,)} = ir J -^5 (w - eg + eg +?± ) (29) 

Intuitively, this follows from the fact that spinons are always excited in pairs and they share the 
momentum transferred from the electron at the interface. At T = 0, the spin liquid is initially in 
its ground state, so we only have contributions from two spin-up spinons that have center of mass 
momentum q±. Equation (29) is the main result of this section, which we shall use to find the 
forms of the spin-current for certain spin-liquids with free-spinon bands in the mean-field picture, 
and then figure out how the spin-current scales with the spin-accumulation voltage V for arbitrary 
spinon dispersions and dimensionality of the system. 
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1. Gapped spinons with quadratic bands 


Let us consider the case of gapped spin liquids in 2 dimensions with a spinon-gap A s , where 
the lowest spinon band has a quadratic dispersion about a minima at k = Q± with an effective 
mass of m*, so that the spinon Green’s function is given by 


G s (k, icJn) — 


lUJr 


- A, - 


(k-Q ±) 2 

2m* 


(30) 


This is true for several ansatz spin liquid ground states 29 ’ 30 , including, for instance, the Q\ = Q 2 
state of the Z 2 spin liquid state on the Kagorne lattice !! ', where the gap and the effective mass are 
given in terms of the mean-field parameters A and Q, and the antiferromagnetic coupling between 
nearest neighbors Jaf by 


A s = 



and 


1 _ 3 J\ F Q 2 


m ’ 


2A, 


(31) 


Equation (29) now leads to the following expression for the structure factor 


S-+ (q±,u) = 


d 2 k 

(2 W 


5 \ u — 2A S — 


(k - Q_ l) 2 (k + q±- Q_ l) 2 


2m* 


2m* 


m 


= —0 w - 2A, - 


*1 

Am* 


(32) 


In general, we may have several spinon bands with minima at different Q± with the same gap A s , 
so we sum over all of them to find the net spin-current via equation (14) in the T —> 0 limit. 


Ispin — 


fJ 2 


J2 fang(k F /Q ± ) J - u)S-+ (q ± , u) 

Q± 


I ) (r - 2A s f 0(V - 2A S ) 


v Qi 


= m (v - 2A S ) 3 @(V - 2A S ) 


(33) 


where we have absorbed all constant pre-factors in rjo to explicitly show the dependence on V. As 
expected, there is a cutoff at twice the spinon gap, i.e, no spin current for V < 2A S , and a power 
law behavior above the threshold. 

Note that in the calculation above, we assume that both spinons come from bands that 
have minima at identical Q±. However, even if they come from different bands, say one with 
minima at Q j_ 1 and the other with < 3 ^ 2 , they will just contribute to add extra pre-factors of 
fang ( &f/(|<2_l,1 + Q. l,2 |) ) in the expression for the spin-current, but would not change either the 
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threshold or the power law behavior. However, if the bands have different spinon-gaps, say A Sj i 
and A <^2 then we expect the spin-current to show a second threshold when the spin-accumulation 
voltage V crosses A Si i + A s ^, as the scattering process then now excite spinons from both bands. 


2. Gapless spinons at Dirac points 


Let us consider spin liquids described by gapless fermionic spinons at discrete Dirac points {Q±} 
in the Brillouin zone. The spinon dispersion is then given in terms of the spinon velocity v at a 
Dirac point at Q± by 


G s (k,iuj n ) — -. (34) 

iu n - v\k - Q j_| 

This is again conjectured to be true for certain spin-liquid ansatz, for example, the 7r-flux state 
of the Heisenberg antiferromagnetic Hamiltonian on a 2d square lattice, which has been argued to 
be stable against U( 1) gauge fluctuations . We again use equation (29) to evaluate the structure 
factor. 


S-+ {q±,oj) = 


d 2 k 
a (2^F 

1 UJ 2 - q \/2 
8nv 2 


6 [uj — v\k\ — v\k + q± 
0 (u; - v\q±\) 


u; 2 -ql 


We now use equation (14) to find the net spin-current for T —> 0. 

_ J 2 A ± u(e F ) 5 _ 5 

hpin ~ non _ 9 ..9 * ~ Vl V 


9607 X 2 V 2 


(35) 


(36) 


where we have again absorbed all pre-factors in r/i to make the H-dependence explicit. The current 
takes non-zero value for any V > 0, as there is no gap to a two-spinon excitation. We have evaluated 
the current for a single Dirac point, although extensions to multiple Dirac points with different 
velocities can be done in an exact analogy with the previous subsection, and will not affect the 
threshold or the exponent in the power law. 


3. Generic spinon dispersions and spatial dimensions 

In this subsection, we are going to generalize the above results for given spinon dispersion in d = 
2 to generic dispersions and arbitrary space dimensions d— 1 of the metal-insulator boundary using 
scaling arguments. Although this approach does not give us the exact-prefactors, it is sufficient to 
find out the characteristic dependence I sp in on V. We would require that the lowest spinon-band 
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has minima at discrete points in the Brillouin zone, which are well-separated from each other. We 
start off with gapless spin liquids with power law dispersions, and find that the exponent of V 
is directly related to the power law in the dispersion and the dimensionality of the system. Our 
results easily generalize to gapped spin liquids. 

Let the spinons have a dispersion given by 

e(k) = v Q \k\ a (37) 

The two-spinon structure factor is proportional to an integral over the allowed phase space consis¬ 
tent with energy conservation. 

S- + (q±,u) ~ J k d ~ 2 dk dfld -2 — v a \k\ a — v a \k + q±\ a ) (38) 

The solutions for k (when the delta function is non-zero) can be written in terms of a dimensionless 


scaling function <&(v a q^_/u) as 

k = q± $(v a q±/uj) (39) 

The delta function in uj can be rewritten as a delta function in k as follows (in terms of another 
dimensionless function <J>i which comes from the Jacobian) 

5{u - v a \k\ a - v a \k + q±\ a ) = 5{k - q _ L <f >(v a q1/u))/ [u Q g" _1 <f>i(u Q g"/w)] (40) 

Now we can see how the dynamic structure factor scales without explicitly evaluating the integral. 

S-+(q±,u) ~ gi~ 1_a ^(v a q1/uj) (41) 

The dimensionless scaling function T must involve a theta function of the form 0(w — C v aQ±), 


where £ is some arbitrary numerical constant that depends upon the exact dispersion. This follows 
from the fact that a large center of mass momentum will inevitably result in a large energy for 
the spinon pair which is precluded by energy conservation. Here, we are assuming that u is small 
enough so that both the spinons come from the bottom of the band(s). 

Finally, we turn to the T —> 0 limit of equation (14) again to find the spin current. 

V r 

(V-u)duj / dq± q 1 ^ 2 d^ d - 2 S- + (q±,u}) (42) 

Because of the © function in S _j_(gj_,a;), the momentum integral is restricted to q < ( u/v a ) 1/,Q: , 

so dimensional analysis tells us that 

J dq± q d ± - 2 dn d _ 2 S-+(q ± ,u;) ~ (io/v a )^ d - 2 - a ^ a 



(43) 
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The integral over oj scales as V 2 , so the final result after putting all this information together is 


I s pin ~ V 2 X y( 2rf -2-“)/“ = V l+2(d-l)/a 


(44) 


As a check, let us see if the scaling matches the previous two exact calculations. In both cases, 
we have d — 1 = 2. For the gapped Z 2 spin liquid in the limit of the gap A s —> 0, we have 
a = 2 and hence, I sp in ~ l/ 1 + 2 ( 3-1 )/ 2 = y3 p or the gapless U( 1) spin liquid with a = 1, we have 


Ispin ~ F 1 ^ 3 " 1 )/ 1 = E 5 . 


For generalizing to gapped spin liquids with a spin gap of A s , all we need to do is make the 
following replacement in all the previous calculations: 


uj —uj — 2 A, 


(45) 


This in turn tells us that the spin current is given by 

Ispin ~ (V — 2A s ) 1+2 ( d - 1 )/° 0(F - 2A a ) (46) 

Equation (46) is the main result of this section. It shows that by measuring the spin current as a 
function of voltage, it is possible to deduce both the nature of the spin gap as well as the effective 
dispersion of the low energy excitations. Note that at the level of low-energy effective field theory, 
the current does not depend on the detailed structure of the lattice, but only on the effective 
continuum dispersion, as expected. 


C. Numerical results for a model Z2 spin liquid state on the Kagome lattice 

In this section, we extend the previous results for a gapped Z 2 spin liquid state via numerical 
calculations. As a model state, we choose the Q± = Q 2 ground state on the Kagome lattice, 
described by Sachdev . The reason for choosing this state for further investigation is that the 
dynamical structure factor measured in neutron scattering experiments on Herbertsmithite single 
crystals is in good qualitative agreement with the calculations in the Q\ = Q 2 ground state by 
Punk et. al . 

Following Sachdev 11 ’, we use a large N expansion technique based on the symplectic group 
Sp(N). To generalize of S~Si to Sp(N), we just extract the part of the Sp(N) invariant scalar 
product S{ ■ Sj that corresponds to In terms of the flavor indices m of the Schwinger 

bosons that make up the spins, it can be written as 

mi ,7712 
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Note that this reduces exactly to S~ of SU( 2) when we have a single flavor. To simplify the 
expression, we note that the N flavors are decoupled in the N = oo mean-field theory, and each 
of the N flavors has an identical Hamiltonian. Therefore, each flavor gives the same contribution, 
which just cancels off the extra factor of N 2 , and we just need to calculate each term for a single 
flavor. The spinon operators that diagonalize the mean held Hamiltonian are linear in the b and 
b t operators, hence the correlation function factorizes as follows 

(S-Sf) = 1 (<6j + <i>!Nt><M.4}) (48) 

Moving to Fourier space and keeping only terms that give contributions to u > 0 after analytic 
continuation, we find that the dynamic susceptibility is given by 

X-+(q±,iUn) = [ u ji(-k)V jm (k + q ± ) + Vji(-k)U jm (k + q ± ) 

X G[( k : TIn, )Grn(k T q±i iuJn T (4:9) 

where q± belongs to the extended Brillouin zone, N s is the total number of sites, U, V are the 
Bogoliubov matrices that diagonalize the mean-held Hamiltonian, and we have implicitly summed 
over all sublattice indices {i,j, l,m}. We are going to use equation (49) to numerically evaluate 
the exact mean-held structure factor. As a side note, we mention that in the low energy limit, 
where k is close to the bottom of a spinon band Q j_, and q± is also small, so that the sum of the 
two spinon energies satishes the energy constraint, we can approximate the elements of the U and 
V by their values at Q j_, and then we recover the dynamic structure factor evaluated in equation 
(32). 

We first plot the momentum-integrated structure factor S _|_(cj) = jr ^^<5_ \~(q,uj) as a func¬ 

tion of energy oj in Fig. 3. We assume mean-held parameters A = 0.695 and Qi = Q2 = 0.4 in the 
units of Jafi which are not self-consistently determined, and lead to a spinon gap of A s ps 0.5. 

We note two specihc features, the jump at w k 0.75 and the peak at w « 1.3. Both these 
features can be understood using the band structure of the spinons for this ground state. The 
spinon spectra has a hat band with = A, and once we have ui > A + A s , we can excite two 
spinons, one of them being at any momentum on the hat band. The second peak presumably 
comes from both spinons coming from the hat band, but is slightly smeared out by the Bogoliubov 
matrices and the hnite width Lorenzian approximation for the delta function in the numerics. If we 
go up to energy scales of V ~ Jaf < ep (this is reasonable as Jaf ~ 200/i for Herbertsmithite , 
but typical €f ~ 10 4 A"), we now can have contributions to the current at large values of Sq x and q±. 


U: i (-k)V* m (k + q±) 
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FIG. 3: Momentum integrated structure factor for the Q\ = Q 2 ground state of the Z 2 spin liquid 

on the Kagome lattice 


In order to investigate the contributions properly, we need to numerically evaluate the spin-current 
starting with the T —> 0 limit equation (6). 

We next plot the spin-current, evaluated numerically, in Fig. 4 as a function of the spin accu¬ 
mulation voltage V. The Fermi liquid parameters chosen for the plot below are kp = 2 (units of 
inverse lattice spacing), and tpj Jaf = 100. 



FIG. 4: Spin-current as a function of spin accumulation voltage for the Q\ = Q 2 ground state of 


the Z 2 spin liquid on the Kagome lattice 
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As expected, we observe the effects of the two features in the dynamic structure factor on the 
spin-current, which is roughly an integral over the structure factor. The step-like jump in the 
structure factor leads to a change in slope in the current around V ~ 0.7, and the spike leads to 
a step-like jump around V ~ 1.3, after which the current saturates. The observation of these two 
distinct features in the spin-current would be strong evidence in favor of the Q i = Q 2 Z 2 spin 
liquid on the Kagome lattice. We note that the Q\ = —Q 2 ground state does not have any flat 
spinon band, and is hence not expected to show any such feature in the spin-current. 


V. CONCLUSION AND OUTLOOK 

In summary, we proposed the use of spin-currents as a gateway to probe the nature of excitations 
in magnetic insulators. Measurement of the spin-current as a function of the spin-accumulation 
voltage can throw light on the dispersion of the low-lying excitations and gap above the ground 
state. In particular, we showed at that in the zero temperature limit, the threshold and scaling of 
spin current with voltages may be used effectively to search for spin liquid ground states in magnetic 
insulators. Finally, we focused on a particular spin liquid ground state, which is a candidate state 
for Herbertsmithite ', and identified some broad features in the spin current which can help to 
identify that state. 

The spin-current is a valuable probe, because once injected into the insulator, the total spin 
is conserved in absence of spin-orbit coupling and random field impurities. We anticipate that 
it may be interesting to study how the presence of disorder in the interface, or the presence of 
non spin-carrying low-lying excitations in the insulator, which couple to the mobile spin-carrying 
modes, (for example, visons coupling to spinons in spin liquids) affect the spin current. 
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Appendix A: Details of spin-current calculations (from II C) 


We begin with the linearized energy transfer u(h, q±, 5q x ) in equations [9], and write the spin 
current from equation [6] as 

7 tJ 2 A ± f d d fa f d d q 

' (27 t)* 


‘-spirit 


4 m 


h\x >0 


Qx kla 


(277) 


n F (%) i 1 - u(n,q±,5q x ))) q x S-+ (q±,u(n,q±,5q x )) 

(Al) 


The integral over \k±\ can now be evaluated, as everything else depends only on the direction n of 
the initial momentum, and the momentum transfer q. Assuming that the density of states v(e F ) 
is approximately a constant close to the Fermi surface, we have 


J >o 
; v{e F ) 


{l-n F {%+V-u{h,q ± ,6q x )) J 


d£i V — u(n, q±, 5q x ) 
Sj 1 - p-P(V-u(h,q ± ,6qx)) 


(A2) 


We can further simplify equation [Al] by getting rid of q x in favor of u. For given q± and n, 
du = Ui?(n • x)d(6q x ) and q x ~ 2 k F (h ■ x), implying 


dg x q x ~ d{5q x )2k F {n • x) 
m m 


—2 du 


(A3) 


This is independent of the direction of initial momentum h. Further, note that the constraint 
q x > kp(h ■ x) is guaranteed to be satisfied by energy conservation, which requires small Sq x . By 

our assumption that S ^(q±,u) is insignificant for large |gj_|, a change of energy due to large 

5q x cannot be offset by another due to large \q±\. The only problem arises when h ■ x is very 
small, but those are insignificant portions of the phase space that we can neglect. Therefore, all 
dependences of the ^-integral on h are removed, and this enables us to do the angular integral. 
Using - >Q we recover the simplified expression stated in equation [10] 


*-spin, 


■kJ 2 A F u{e F ) 


du d d 1 q 1 _ (V - u) c ^ ^ 

2ir (2TT) d ~ 1 1 - + 


(A4) 


The calculation for I spin,i (equation [11]) is analogous, with the only change coming from the 
different occupancies of the initial and final states. 

We now discuss the case when the dynamic structure factor has a minima at large transverse 
momentum Q±. The trick is to note that although the momentum transfer can be large, the 
energy transfer at low temperatures is always small, i.e, u < V <C e F . Therefore, we can expand 


in small parameters about the point of elastic scattering. To do so, we first solve for a longitudinal 
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momentum transfer q x o which satisfies where Q = q x qx + Q±. 

% - %-Q = 2kF id -L • fl + 9x0(« • 5)) - (<2i + = 0 

1 /2 

=> feo = (n • x) + (k 2 F (h ■ x) 2 + 2k F (Q± ■ h) - Q 2 ^ (A5) 

The contraint q x q > k F (h • x), required for reflection, implies that only the positive square root 
can contribute. As q x is real, only some values of h are relevant. Specifically, we require 

k F (h ■ x) 2 + 2k F (Q± ■ n) > Q\ (A6) 

We need to evaluate the angular integral over angular regions of the Fermi surface consistent 
with the above constraint. If Q±, the magnitude of the transverse scattering wave vector, is too 
large compared to the Fermi momentum k F , then there is no scattering consistent with energy 
conservation, and hence there is no spin current due to this process. 

Now, we can expand about the solution for elastic scattering for small a;, and keep only linear 
terms in 8q x and q±. 


q = ((ho ~ 8q x )x + Q±- q± 
1 


u = — 
m 


_i_ - k F h ■ q _l 


+ 0(5ql,ql) 


(q x o - k F (h ■ x)) Sq x + (Q 
where we used 2 k F (Qx ■ n + q x o(h ■ x)j — ( Q]_ + q x0 ) = 0 


(A7) 


We revert to our previous formalism, and replace the integral over q x by an integral over a 
with the only change being in the pre factor appearing the angular integral. For fixed q± and h, 
duj = ^ (q x o - k F (h ■ x)) 5q x , and q x ps q x0 , implying that 


dq x q x _ d(5q x )q x o _ ^ / QxO 
m e m e \k F (h-x)-q x 0 


1 + 


k F (h ■ x) 


y [k 2 F (n ■ x) 2 + 2k F (Q±_ ■ n) - Q 2 ± ^j 


\ 1/2 


duj 


Note that for \Q±\ <C k Fl we get back our previous result which corresponds to scattering at 
Q l = 0. This acts as a check on the above calculation, and also shows that the calculation can 
be generalized as long as we have low energy excitations in the spin-system about a set of isolated 
points in momentum space which are well-separated in the Brillouin zone. 

The factor multiplying du will change the result of the angular integral over initial momenta, 
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but the remaining calculation remains unchanged, and we have 

/ 


tQ± = A±v(e F ) 

1 spin,']' 


dfl 


^ Jk^(n’x) 2 +2kjr(Q±-fi)>Q‘^ $d— 1 

dto d d ~ 1 qj_ ( V-u) c N 

27 t {2n) d ~ 1 1 - e-d(y-^) + [Q±,UJ) 


1 + 


k F {n ■ x) 


y f k 2 F (h ■ x) 2 + 2k F {Q± ■ n) - Q 2 


1/2 


TrJ 2 A ± u(e F ) , . f du d d 1 q± {V-u) 

-i- /«w(*f/Gi) J 2tt ( 27r) rf ~i 1 - S- + (l±M 


(AS) 


where f a ng(k F /Q± ) is the angular integral referred to in equation [13]. Typically, k F and Q± have 
the same order of magnitude, and then the angular integral is an overall factor of 0(1) (the exact 
value is determined by the constraints set by the ordering wave vector Q±). 

Taking into account that there can be multiple such minima in the dynamic structure factor 
at large finite momenta and scattering to momenta patches around these minima are inde¬ 

pendent as long as the minima are well-separated, we arrive at equation [12], stated below for the 
sake of completeness. 


I, 


spin,']' 


irJ 2 A±i'(e F ) 


£/■ 
Q 


ang 


(k F /Q±) 


du d d l q_ l (V — u) 

2n~ 


(2ti tY 1 - 1 1 - e ~P( v ~u) S ~ + 


(A9) 


The expression for I S pi n ,i follows in exact analogy to the above calculation. 


Appendix B: S _ (q±,u) for an antiferromagnetic interface (from III A 2) 

We evaluate the dynamic structure factor for a Neel-ordered state on a d-dimensional cubic 
lattice using the Holstein-Primakoff transformation. First, let us consider the case when the Neel 
vector points parallel to the spin-quantization axis in the metal (chosen to be z). We have up-spins 
on sub lattice A and down-spins on sub lattice B, with total number of spins being be 2 N, and the 
coordination number of each spin is z = 2d. Therefore we define 

i G A, Sr = 4(25 - ofa*) 1 / 2 ; 5+ = (25 - a}ai) 1/2 <H, and 5? = 5 - ojo* 
i € B , 5+ = 6j(25 - b\bi) 1/2 ; 5+ = (25 - bt&i) 1/2 k, and 5f = -5 + b\k (Bl) 

and do an expansion in 1/5. The Heisenberg Hamiltonian Haf = Jaf can be written 

in terms of the Holstein Primakoff bosons as 

Haf = -JAFNS 2 z+J A FSzY^W % a % +b\b % + 1% {a$_ % +a\b i _ J .)\+O{S 0 ), where 7jg = * ^ e lkS 

fc S^n.n 

(B2) 
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This can be diagonalized by a Bogoliubov transformation, using 

a k = u k a k+ v k^U ’ b k= u kPk + v k a U 

with Ur = u_r = cosher), v £ = v_r = sinh(0r) and tanh(2dr) = —'yr (B3) 

The Hamiltonian is diagonal in terms of the Bogoliubov quasiparticles, 

Haf = -JafNS 2 z - JafNSz + E^ala^ + /?!/% + 1), E % = JafSz^Ji - 7 S 

fc 

S _| .(q±,uj) may now be calculated from definition using the expression for the spin operators in 

terms of in terms of the quasiparticle creation and annihilation operators. After some algebra, we 
find 

S-+(q±,u) = 2tt S(u$ ± +v $ ± ) 2 [d(w - E$ ± )(l + n(P$ ± )) + S(u + Ef ± )n(a-$ ± )\ (B4) 

At T = 0, only the delta function with positive u> contributes to the spin current as there are on 
quasiparticles initially in the system. For low momenta, we have 

Eg ± ~ v s \q±\, and (u$ ± + v^) 2 = cosh(2f% x ) + sinh(2% x ) = , j 7<U = ~^= (B5) 

\ 1 r 2yd 

which leads to the following expression for the dynamic structure factor for the T = 0 antiferro- 
rnagnet 

S-+(q±,v) = “ v sQ±) (B6) 

For the case when the Neel order is perpendicular to the spin quantization axis in the metal, we 
assume that spins on sub lattice A are pointing in the y direction and the spins on sub lattice B are 
pointing in the — y direction. In this case, we can still use the Holstein-Primakoff representation of 
spins after doing a 7r/2 rotation of our coordinate system with respect to the x axis. In the rotated 
coordinate system XYZ we have X = x,Y = —z, and Z = y. Remembering that our original 
definitions of were with respect to the old axes, let us denote our spin operators by S in the 
new set of axes. Then 

S ± = S x ± iS y = Y x ± iY z (B7) 

We can now express these in terms of the usual Holstein Primakoff bosons, and after some algebra, 
find the following dynamic structure factor in the large-S approximation 

2tj- S 

S-+i.q±,u) = -^{Uq^+Vqj 2 [<S(w - E ?± )( 2 + n(a ?± ) + ra(%)) + <5(w + E ?± )(n(a_ ?± ) + ri(/3_ fx ))] 

(B8) 
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Using the low momentum limit from [B5] and taking T —> 0, we arrive at the expression in equation 
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